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ABSTRACT 


Based upon simultaneous observations of neutral-lapse conditions with 
the corresponding surface wind value, both for ship Victor, values of 
the roughness parameter have been computed. For these computations 
an adjunct relationship for drag coefficient after Kung (1963), is 
also employed, the latter relationship being primarily statistical 

in nature. It was found that the roughness length undergoes a 
variation of ten orders of magnitude as the ratio of the observed 
Surface to geostrophic winds varies through its complete sample 


range from approximately 0.48 to 0.95. 
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SYMBOL 
Us, 

H 

L 

V 

V20 

Vg 

Z 

ZO 


vm 
Cp 


ETS0F sYMBOLS 
MEANING 
friction velocity 
turbulent heat flux 
scale length (see introduction) 
wind speed 
wind speed at 20 meter height 
surface (1000 mb) geostrophic wind speed 
height above surface 
Surface roughness parameter 
averaged stress vector exerted at sea-air interface 
geostrophic drag coefficient 
air density 
von Karman constant taken as 0.4 
dry bulb temperature 
sea surface temperature 
Surface Rossby number 
coriolis parameter 
potential temperature 
Specific heat of air at constant pressure 
acceleration of gravity 
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1. Introduction. 

According to the accepted theories on atmospheric turbulence, 
including those dealing with the Monin-Obukhov similarity principles, 
there exists a turbulent, atmospheric boundary layer characterized 
by certain near-constant parameters. These parameters are Uy, the 
friction-velocity, the turbulent heat-flux H, and a scale-length 


L defined by 
3 


U 
L=-7ga (1) 


Oo Of 
In the present study, where case selection was restricted to neutral- 


lapse cases at Ship Victor, L is infinite (H=0) and the usual logarithmic 


profile 


pal, = 


has been assumed to be valid. Here V is the wind speed at the anemometer- 
level, which was taken as Z = 2000 cm (after some correspondence with 
cognizant U.S. Coast Guard Officials). No attempt was made to introduce 
a second statistical degree of freedom, in the form of the datum-level 
displacement, mainly because accurate wave-heights were not included in 
the coded data cards for any of the available periods under study. 

As is well-known, most operating numerical prediction centers, for 
example, that of the National Meteorological Center in Suitland, Maryland, 
make use of average values of the geostrophic drag coefficient, Cp, 


obtainable from 


T = oe VY," = ¢ Ue (3) 


The right side of (3) is dependent upon the roughness parameter Zo 
(see equation (4) below), grossly averaged, and then taken to be 
representative of NMC grid points. Operationally, there seems to be 
little cure for injecting daily deviations of Cy from its seasonal 
average, but this study was undertaken to investigate the possibility 


of the existence of such deviations. 


2. Selection of the Data. 

All data were taken from punched data cards provided by the 
National Weather Records Center at Asheville, North Carolina. The 
data cards employed spanned the period 1959 - 1960 at Ocean Station 
Vessel "Victor" employing only the section 116 (Marine Surface Obser- 
vations) part of the coded data. Since only 6 to 10 meters of elevation 
is involved between sea and instrument levels, a neutral lapse would be 
"isothermal" over this vertical span. However, since both the sea and 
dry bulb temperatures have, with very few exceptions, been preliminarily 
rounded off to the nearest degree, the selection criterion for "neutral" 


was relaxed to the following: 
7 = IE = 1+t/'F 


Each of the neutral cases were carefully cited as to date and time, 
and if the observation time coincided with that of a 00 GCT sea-level 
analysis in the file of Northern Hemisphere Synoptic Weather Maps 
covering the same period (1959 - 1960), the particular neutral case was 
accepted for further computation. 

For each case accepted, the shipboard wind speed in knots was 
recorded from the card file. For the same case the geostrophic wind was 


carefully read off by measuring the perpendicular spacing between 5 mb 


isobars on the synoptically=timed sea-level chart. These measurements — 
were made using dividers; then appropriate distance measurements were 
converted into a geostrophic wind at the latitude of "Victor" (34°N, 
164°E). It was found that use of the geostrophic wind scale on the 
individual weather maps yielded objectively sound values provided 

cases where (1) fronts occurred close to Victor, or (2) the isobars 

in the vicinity of Victor were strongly-curved cyclonically, were 


rejected. 


3. Derivation of the roughness-parameter equation. 
First consideration was given to the neutral wind profile applicable 
at 20 meters: 
Vi oe 2000) (4) 
k x, 
with Ux = 7 $F considered constant in the surface layer. Kung (1963) 


defines Ux of equation (4) in accordance with 


Ux. we yy Vg (5) 


where Cy is his drag coefficient, but is to be interpreted as the 
Square root of that of Cressman (1960). Furthermore, Kung's statistical 
analysis of Ux resulting from (4) has led to the following empirical 


result for Cp: 


— _ v 
- Leq Ko, -—O.556 ya Fz, ye 


It should be noted that Ro, as defined in equation (6) is a dimensionless 
parameter describing the large-scale frictional aspects of the boundary 


layer, and has been called the surface Rossby number. 


Combination of equations (4, 5) leads to the result: 


Vao— 


G= ACY) _ k (7) 
Lr (=22") 2.3026 Day (222 ) 


This result, combined with equation (6) yields 


o. 205 | Lvs 
A ~ lay Z ~hy F- OSS6 (2.3026) [hy 2000 - Leg Z/ 


It should be noted that (8)is valid only for the statistical range of 





Kung's resuit (6) and for those conditions under which a neutral 
boundary layer may be considered to apply. Within these rather broad 


limitations, the solution for log #o is 


: - a yO *) es a 4 Bt “ A, , ) 
3 sa = Ly 4 a F-OSS6 — os ve - . din acter = or ) (9) 
< . 1H) Q 


or more specifically, with k = 0.4, log 2000 = 3.30103 and 


log f = - 4.0886, equation (9) becomes 


: Yo AD LNG (10) 
z= Leg Vg + 3,536 = Con A 0.2 516) 
Lag 0 y (/- OP Y Fy 44, ) 


Introduce the notation 


x = SS (11) 
(1- 0.9474 VXy) 


and the resulting expression for log #0 finally becomes 


ey 8, = Ay \y 4 2.5226 - Gay Vy 40.2316) X (12) 


Since it is customary to solve for Zo in cm, and f has been expressed 
in sec “4, all wind speeds V20> Vg have been converted from kts to cm 
ial using the conversion factor 

1 knot = 51.479 cm sec”? 


Equation (12) is the working relationship for the calculations made in 


this study. 


4. Results. 

All pertinent boundary layer parameters have listed in Table 1. 
(see Appendix) for each of the 55 cases which fitted the acceptance 
criteria. It is noted that Va9/Vg varied from a minimum of 0.489 on 
1 January, 1959 to a maximum of 0.969 on 25 November 1959. Moreover, 
Vo0/Vg appears aS a linear variable within the denominator of X. It 
is therefore evident that the single factor which most strongly affects 
the value of log Zo is X, which varies from a minimum of 1.708 to a 
maximum of 5.596, corresponding respectively to the two wind ratios 
just cited. Qn the other hand, the contribution of the only other 


variable in equation (12), log V., was relatively minimal varying 


g 

only by a multiplicative of factor of 0.90 to 1.10 relative to the 

sample-mean of log Vg: 
Thus each value of log #0 was coupled with the simultaneous ratio 

of X, defined by equation (11). As the above comments suggest, there 

was an extremely regular tendency for log Zo to decrease with increasing 

(V29/Vg). This was true both in the 55 individual computations, as 

well as in class=grouped averages of both log Zo and Vag/ Vg: The 

grouping was based on classes of V0!“ according to the values of these 

wind ratios in the number ranges 0.45-0.50, 0.50-0.55,......... ; Uroo-0.290, 

0.90-0.95. The corresponding class—mean averages of log Zo were then 


computed for each wind-ratio class. The paired class-average values 


were as follows: 





+0.4500 {+0.1986 j-0.8296 |-1.3031 |-2.2320 










0.75-0.80] 0.80-0.85]0.85-0.90}0.90-0.95 
0.7651 0.8202 0.8624 0 9ieil 
-3.4514 |-4.5050 |{-5.7238 |-8.6549 


The results by class-groups so closely reflected the individual 


(log Zo 
Average 


pairings that the final results have been graphed in Figure 1, solely by 
paired class-averages. The resulting graph essentially presents log 2Z0 


in the functional form 


fag, = (GV, +3,5326)-[Gyi, ¢ 0.2316] (/- o.00 74 Ma) (13) 


where log Mig is the sample-mean of the individual values of log \y, In 


other words, log Zo has a general simplified form 


Z =/ 
bay B= 0, - Cx (1-094 79) (14) 


where C, is the constant within the first parentheses of (13) and Ca 
that within the bracket of equation (13). 

Figure 1 shows that log Zo tends to decrease towards=@@ as V50 
approaches Vg whereas at the other extreme for Vog/Vq = 0.489, #0 = 5.5 cm, 
an unusually high value for an oceanic area. Upon exclusion of that 
class-group Vo9/Vg with values in the range 0.95 to 1.00 (these were 
considered unrealistically large for the neutral, steady state case), 
the mean-value of Zo turned out to be intermediate between the extreme 
oceanic roughness-values Zo = 0.1 cm (for V2. 10 kts) and Zo = 0.01, 
for V.< 10 kts, employed by the General Circulation Laboratory of ESSA. 


Their approach, also empirical, evidently presumes that stronger surface 


winds and a consequent greater roughness parameter, ensues from the 


assumption of a greater geostrophic wind. Such a modeling assumption 
need not be inconsistent with the results found here, but merely an 
expression of the relative scopes from which the turbulence problem 
has been viewed. Some further physical justification of equation (14) 


is offered in Section 5. 


5. Conclusions 

It has already been shown in the preceding section that the 
application of neutral boundary layer wind data chosen purely at random 
leads to a tenfold variation in Log #0. However, recent theories of 
the planetary boundary-layer, specifically those of H.H. Lettau (1961) 
and Blackadar (1962) for the barotropically neutral atmosphere, but 
with variable exchange coefficient K=K(Z) led to spiral-type solutions 
with increasing 2. One of the main shortcomings of these theories 
in ascribing proper values of the surface-stress is the assumption of 
zero acceleration in the large scale wind field. Kung (1963) has given 
average values of V16/Vq at numerous geographic locations over the 
oceans, these based upon rather long-averaging periods. When one averages 
over his summer and winter climatological ratios, the resultant V16/Vq 
reflects both near-neutral and quasi-steady state conditions. This average 
value of Vi6/Vg was 0.68, and extrapolation gives V20/" = 0.70. 

In this study, our median class centered about 0.70, but large 
deviations in V20/Vq relative to 0.7 did occur, with higher wind ratios 
reflecting supposedly a smoother air-ocean interface. The interpretation 
of this surprising result must lie in the use of individual data applied 
to steady-state theory, which represents improper application of limited 


theory. 


The defects seem to occur in the areas of 

(1) the non-inclusion of acceleration terms in the boundary-layer 
theory. 

(2) the slow variation rate of a given sea state, even though the 
impressed large-scale geostrophic wind undergoes significant 
change. 

Both areas (1) and (2) are interrelated through the basic unknown, 

the large-scale acceleration, and its feedback rate to the interface. 
For example, consider our Vo9/Vg = 0.5 relative to steady-state 


turbulence theories. This was taken to mean that V, (at gradient level) 


aul 
had responded to an acceleration to which the surface wind was unable 

to respond, except at much greater lag time. Hence in entering equation 
(12), X is below average and the resulting log #o greater than average. 
On the other hand, an individual case wherein V29/Vg = 0.90 compared to 
the steady-state value of 0.7 indicates that overhead a decrease in Vg 
had already taken place without a corresponding decrease in Vog> due 
largely to the greater lag-response time at the surface. 

The main conclusion of this study is not that of questioning the 
validity of values of #0 in operational use. Rather, its purpose is 
to point to areas where the turbulence theory seems to be in need of 
improvement, especially with regard to the inclusion of unsteady-state 
conditions in the free-atmosphere, and the difference in response-time 
in the boundary layer as the interface is approached. The lag time 


consideration is especially compounded when there is an oceanic 


interface subject to storm-wave disturbances of short wavelength. 
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APPENDIX I - TABLE 1 


Table of Variables Needed In Solving 
For Log #0, Listed by Case and Date 
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Date 


2) 39 
13/59 
ye 759 
ofl 59 
23/59 
25/09 
aby O92 
S07 59 
apc /59 
2/3/59 
2/8/59 
2/10/59 
Cia 09 
aye / 59 
2/13/59 
@73/09 
By a / 59 
Sy 59 
Oy Es a, 
S731 59 
4/4/59 
4/5/59 
4/9/59 
4/10/59 
4/11/59 
4/14/59 
4/17/59 
4/18/59 
ey 37 59 
10/18/59 
10/20/59 
10/28/59 
0729/58 
10/30/59 
ny? £59 
yaa 59 
1 / 59 
22759 
ly 25/59 
yey 739 
12/14/59 
27/157 59 
12/16/59 
12/17/59 
12/18/59 
12/20/59 


T.(°F) T(OF) Vogt kts lap, Log V, 
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63 
67 
6] 
6] 
60 
60 
60 
98 
oo 
oe, 
62 
6] 
63 
63 
64 
63 
64 
63 
62 
63 
aS) 
6] 
63 
64 
64 
63 
63 
63 
74 
ic 
WZ 
65 
69 
69 
65 
7] 
71 
6/ 
69 
70 
68 
68 
68 
68 
69 
69 


64 
67 
60 
6] 
6] 
6] 
60 
58 
98 
60 
62 
60 
63 
63 
65 
63 
63 
62 
61 
63 
og 
6 | 
63 
64 
64 
64 
63 
63 
74 
7| 
is 
64 
68 
68 
65 
7) 
71 
66 
69 
70 
67 
68 
69 
68 
69 
69 


65 
90 
90 
23 
47 
ZO 
90 
28 
f(a 
4] 
Zo 
E3 
Md 
14 
i] 
43 
30 
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23 
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3] 
3] 
29 
e 
38 
38 
3] 
a 
2/ 
24 
40 
29 
Jil 
26 
17 
aT 
34 
33 
38 
Si 


TABLE 1 
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22 sm 73 
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| 4 (ohehs) 
20 gmc03 
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26 ohn thks, 
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is S073 
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ZS Ripe) 
2 | Salen 
30 30S 
2 | 32208 
2° 3.174 
“| bie 40/8: 
25 ein 4S)! 
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Zi 32205 
(as 3.174 
18 3 ale 
Zz Saag 
32 3.314 
28 3.4 
30 31203 
24 Selig. 
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Zc 3.143 
Ze 3.243 
26 Sao 
656, o)5 25. | 
Si gecol 


Wi 


Vo0/Vq 
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le, 
56] 
Jeo 
seyey 
. 863 
steken/ 
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soe! 
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.677 
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pres | |e 
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Boy) 
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ale 
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soul 
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SENS) 
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1 Og 
448° 


Log #0 


le 
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+0 
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Boe 


-0 


ay. 
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Fig.1 Graph of Logzo against Vao/Vg for Ocean Vessel Victor 
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